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Abstract 

Gravitational field mismodeling produces errors in the estimated orbital 
motion of near Earth satellites. This effect is studied using a linear perturba- 
tion approach following the analysis of Kaula. The perturbations in the orbital 
position as defined by either orbital elements or Cartesian components are de- 
termined. From these perturbations it is possible to ascertain the expected 
signal due to gravitational mismodeling that would be present in station-to- 
satellite laser ranging measurements. This expected signal has been estimated 
for the case of the Lageos satellite and using the predicted uncertainties of the 
GEM-T1 and GEM-T2 gravity field models. The results indicate that observ- 
able signal still exists in the laser range residuals given the current accuracy of 
the range measurements and the accuracy of the gravity field models. 


(NASA-CR-196057) SATELLITE LASER 
RANGING AND GRAVITY FIELD MODELING 
ACCURACY Final Report (Colorado 
Uni v. ) 72 o 


N95-12837 


Unci as 


G3/13 0026407 


’Department of Aerospace Engineering Sciences, Campus Box 429, University of Colorado, Boul- 
der, Colorado 80309 



1 Introduction 


Satellites in Earth orbit are affected by a wide range of gravitational and non- 
gravitational perturbations. An accurate modeling of these perturbations is required 
if the full scientific benefits of geodetic satellites, such as Lageos, are to be obtained. 
This modeling is sometimes limited by a lack of understanding of the physical pro- 
cesses producing the perturbations, but it is more often limited by the accuracy of 
the parameters that are needed in the model definition. A particular example of this 
latter scenario is the modeling of the Earth’s gravity field. The mathematical form of 
the gravitational perturbation due to the Earth is well defined (using spherical har- 
monics) but the constants used in the model are of varying accuracy. This leads to a 
mismodeling of the gravitational effect and a resulting mismodeling of the predicted 
satellite motion. This mismodeling is observable if accurate tracking measurements 
of the satellite exist; such as laser range tracking. 

To ascertain the level of mismodeling of the gravity field existing in current models, 
a linear perturbation analysis has been applied. This approach makes it possible to 
predict the expected magnitude of orbit error due to gravity field modeling error and 
it also has been extended to determine the expected signal that would be observed 
with laser range tracking. 

The application of using spherical harmonics for modeling the Earth’s gravity 
field is first summarized. This is followed by a discussion of satellite dynamics and 
the resulting orbital motion due to the presence of a nonspherical gravity field. The 
perturbation approach, based on the method of Kaula, is then presented. Pertur- 
bations in orbital elements and satellite position are included. From these results, 
the expected perturbation in a station-to-satellite range measurement is then derived. 
These linear perturbation results facilitate a covariance analysis based on the expected 
uncertainties of present day gravity field models. Numerical results are presented for 
the case of the Lageos satellite and the GEM-T1 and GEM-T2 gravity field models. 




2 Gravitational Field 

2.1 Spherical Harmonic Model 

The common approach for modeling the gravitational field of a planetary body is 
through the spherical harmonic representation, 

U = -J2 (— ) £ p tm (sin <t>)[Cim cos mX + Stm sin mA] (1) 

r 1=0 ' r ' m *0 

where: /z is the product of the universal constant of gravitation G and the mass of 
the Earth M; a e is the semi-major axis of the Earth’s reference ellipsoid; r, 4 > , A are 
the satellite distance, latitude, and longitude, respectively, in a body-fixed coordinate 
system; C* ro , S< m are spherical harmonic coefficients of degree i and order m; and 
P tm are the Associated Legendre Functions of degree l and order m. A gravitational 
model consists of a set of constants that specify n, a e and the C( m , S( m coefficients. It 
should also be noted that such a set of constants also implicitly defines a body-fixed 
coordinate system. The coordinate system defined is precisely that which was used 
in the solution of the spherical harmonic coefficients. 

This representation of the geopotential can be thought of as consisting of three 
constituent parts, 

U = U 0 + Ui+U 2 (2) 

The first part is simply the leading term of the expansion corresponding to the degree 
and order zero term. The Associated Legendre Function, Poo has a value of one as 
does the Coo coefficient. So the leading term is simply, 

U 0 = - (3) 

r 

This is the familiar potential resulting from treating the body as point mass and that 
used for deriving the fundamental results of two body motion. 

The second part of the spherical harmonic representation are those terms (besides 
the above two body term) which do not have a longitude dependence. These are 
the terms corresponding to m = 0 and are denoted as the zonal contribution to the 
potential, 

= ( 4 ) 

The degree 2 zonal term models the contribution due to the planetary oblateness. 
As such, it is the second largest contributor to the overall potential following the 
central body contribution. (The degree 1 term is zero assuming that the center of the 
Earth- fixed coordinate system coincides with the center of mass of the Earth.) The 
notation Jt is often used for the zonal coefficients instead of the above Ct,o- The two 
notations simply differ in sign, 


J( = —C(q 


( 5 ) 



( 6 ) 


so that the zonal part of the potential could also be written in the form, 

r s 

The C notation will be used throughout this report. 

The remaining part of the spherical harmonic representation is that part depend- 
ing on longitude, 

U 2 = - (— ) H P(m ( sin <t> ) [Ce m cosmA + S tm sin mA] (7) 

r £A r ' 

The largest longitudinal contributor to the potential is usually the degree 2 and order 
2 terms. These terms represent the amount that the planet is “out of round” about the 
equator. (As with the degree 1 zonal coefficient, the degree 1 and order 1 coefficients 
will be zero under the assumption that the center of the coordinate system coincides 
with the center of mass.) 

The spherical harmonic representation of the potential (equation (1)) can then 
also be written as, 



In the general case, where temporal variations in the potential exist (e.g., tides), 
the spherical harmonic representation is still valid though the geopotential coefficients 
( C tm , Stm) then become time dependent. 

2.2 Normalized Spherical Harmonic Model 

The spherical harmonic coefficients ( C tm , S (m ) appearing in equation (1) are unnor- 
malized. These coefficients tend to very small values as the degree increases. This 
is partly a consequence of the nature of the Earth’s gravity field but is for the most 
part due to the fact that the Associated Legendre Functions tend to large values as 
degree increases. Thus it is numerically advantageous to normalize the Associated 
Legendre Functions and the coefficients. The normalization is achieved by multiply- 
ing the Legendre functions by a scale factor depending on the degree and order of 
the function. The resulting normalized Associated Legendre Functions maintain the 
same magnitude with increasing degree. Denoting normalized values by an overbar, 
the normalized Associated Legendre Functions are, 



P (m = (2-<U)(2£+l) |^™jj 


(9) 



where the Kronecker delta, 5 m o, is equal to 1 if m is zero and equal to 0 if m is greater 
than zero. The geopotential coefficients, Ct m and 5< m , are normalized by the inverse 
of this scale factor, 


C tm 1 = 1 (l + m)! l 1/2 f Ctm \ 

j L( 2 “W( W+ l)(^“ m )!J \ Sim J 


( 10 ) 


The spherical harmonic expansion of the geopotential (equation (1)) can now be 
written in terms of normalized quantities, 


OO / v f t 

u = ( — ) £ Am (sin 0) [^ m cosrnA + 5 <m sinmA] (11) 

r 1=0 ' r ' m =0 

This is usually the preferred formulation for numerical implementations of the spher- 
ical harmonic representation. For some analytical developments it is simpler to work 
with the unnormalized form (equation (1)). 


2.3 Earth Gravity Field Model 

Numerous spherical harmonic models of the Earth’s gravity field have been developed. 
These models are primarily based on the Earth-based tracking of low-Earth orbit 
satellites. Other data types that are valuable in estimating the Earth’s gravity field 
include surface gravity measurements, satellite-to-satellite tracking and more recently, 
satellite radar-altimeter measurements of the ocean surface. The maximum degree 
(£) of the spherical harmonic representations of the Earth is more than 300 in some 
models. Models based solely on satellite tracking data usually have a maximum degree 
of approximately 50. Indeed, for most satellite applications, high degree models are 
not needed due to the insensitivity of the satellite motions to the small scale features 
represented in such models. 

A good general purpose model for satellite applications is the GEM-Tl model 
developed at the NASA Goddard Space Flight Center (GSFC) [Marsh et al., 1988]. 
This model is based solely on the tracking data of Earth satellites and is complete 
to degree and order 36. The accuracy of this model is greatly improved with respect 
to earlier Goddard Earth Models (GEM). In addition to estimating the spherical 
harmonic coefficients of the gravity field, GSFC also estimates the accuracy of those 
coefficients. Such accuracy estimates are very valuable when attempting to estimate 
the orbit error which may be induced when using the model for orbit propagation. 

The GEM-Tl model was followed in development by the GEM-T2 model [Marsh 
et al., 1990]. This model was a result of continual refinement in analysis techniques 
as well as the inclusion of more and newer observations. The GEM-T2 model is 
similar to the GEM-Tl model in that it is based solely on the analysis of satellite 
tracking data. The GEM-Tl and GEM-T2 models were the primary models used in 
the analyses presented here. 


Table 1. The GEM-T1 Gravity Model^ 
Normalized Zonals Cm in units of 1 x 10~ 9 
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The GEM-T1 normalized zonal coefficients, complete to degree 36, are given in 
Table 1. The nonzonal coefficients are listed in Table 2. Along with the coefficient 
values in each table, the estimated uncertainty (<r) of the individual coefficients is 
also given. The gravitational constant and equatorial radius specified for the GEM- 
Tl model are, 

H = GM = 3.98600436 x 10 14 m 3 /s 2 and a e = 6378137. m (12) 

Several points can be made by examining the coefficient values and their uncer- 
tainties. Foremost, the value of the second degree zonal coefficient is seen to be more 
than two orders of magnitude larger than any other coefficient. The next largest coef- 
ficients are those of degree 2 and order 2. Analogous to the second degree zonal which 
represents the oblateness of the Earth, these coefficients correspond to the ellipticity 
about the equator. 

Also evident is that the magnitude of the coefficients decreases significantly as the 
degree increases (keep in mind that these are normalized coefficients and effectively 
have equal weight in their total contribution to the gravitational potential). This 
characteristic has been formalized in the so-called “Kaula’s Rule” \Kaula 1966]. This 
rule gives the expected size of the Earth's normalized harmonic coefficients of degree 
t to be ±10 ~ S /P. This rule of thumb allows one to estimate the expected magnitude 
of a gravitational coefficient if a value is not otherwise known (this is particularly 
convenient for higher degree coefficients for which accurate estimates have generally 
not been obtained). The magnitude of the GEM-T1 coefficients (complete to degree 
and order 36) are plotted in Figure 1 along with Kaula’s Rule. 

The estimated accuracy of the various coefficients shown in Tables 1 and 2 show 
that the lower degree coefficients are the best determined and the accuracy degrades as 
the degree increases. This variation in accuracy is a reflection of the fact that satellite 
tracking data was used to solve for the coefficients. As discussed in later sections, 
the sensitivity of the satellites to the harmonic coefficients decreases as the degree 
increases. That is, the low degree coefficients produce large perturbations to the 
orbital motion and the high degree coefficients produce much smaller perturbations. 
The ability to recover high degree coefficients is a direct function of the accuracy of 








Figure 1: Absolute value of the normalized GEM-Tl coefficients and Kaula’s Rule. 

the tracking data being utilized and also the geographic distribution of that data. 
Since the high degree coefficients represent fine scale features in the gravity field it is 
necessary to have wide geographic coverage, and accurate tracking data, to completely 
capture such details. The low degree coefficients on the other hand represent the large 
scale features (continental in size) and it is possible to accurately model such details 
given sparser geographic coverage. 

Overall, the accuracy of the geopotential acceleration is dependent on the con- 
stants Cf m and Stm and the accuracy of the coordinate transformation from body- 
fixed to inertial coordinates. The accuracy is also clearly dependent on the degree of 
truncation of the infinite series describing the potential. It should be noted that the 
quantity a e enters the potential strictly as a scaling factor and thus does not affect 
the accuracy of the geopotential computation. Of these possible error sources, the 
accuracy of the C im and S tm coefficients is currently the limiting factor in precise 
low-Earth orbit determination. 



T^ble 2. The GEM-T1 Gravity Model. Normalized Sectorials and Tesserals 
Ct™ and S/r„ in units of 1 x 10~ 9 
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Table 2. Continued 
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11 

6 

8.47 

5.0 

24.29 

5.1 

12 

5 

30.60 

8.0 

- 1.47 

8.1 

12 

6 

1.39 

6.2 

45.83 

6.1 

13 

5 

59.62 

10.0 

57.45 

10.5 

13 

6 

— 22.39 

7.7 

- 11.84 

7.7 

14 

5 

22.80 

10.1 

- 11.61 

10.1 

14 

6 

- 3.19 

9.2 

6.51 

9.1 

15 

5 

16.07 

11.5 

10.89 

12.1 

15 

6 

27.23 

10.5 

- 51.71 

10.4 

16 

5 

- 7.72 

11.3 

- 1.68 

10.9 

16 

6 

17.99 

11,1 

- 26.78 

11.0 

17 

5 

- 11.15 

10.6 

- 5.60 

11.2 

17 

6 

.29 

12.2 

- 30.42 

12.0 

18 

5 

1.74 

11.7 

21.13 

11.4 

18 

6 

31.20 

11.0 

- 8.56 

11.0 

19 

5 

- 2.45 

10.7 

17.38 

10.8 

19 

6 

- 6.25 

11.1 

3.92 

10.8 

20 

5 

- 10.42 

10.9 

.30 

10.7 

20 

6 

12.76 

10.6 

.97 

10.6 

21 

5 

17.76 

11.1 

- 15.89 

11.2 

21 

6 

4.21 

10.1 

- 8.35 

10.0 

22 

5 

- 4.63 

10.5 

-.13 

10.5 

22 

6 

14.63 

9.7 

2.43 

9.8 

23 

5 

1.99 

10.8 

- 7.94 

10.8 

23 

6 

10.00 

9.7 

4.96 

9.5 

24 

5 

- 14.08 

9.7 

- 7.98 

9.7 

24 

6 

-.31 

9.5 

-.66 

9.6 

25 

5 

- 2.46 

9.5 

- 2.35 

9.6 

25 

6 

5.94 

9.6 

- 6.74 

9.5 

26 

5 

4.30 

8.9 

10.65 

8.9 

26 

6 

8.54 

9.0 

3.20 

9.1 

27 

5 

- 1.59 

8.7 

3.71 

8.7 

27 

6 

1.77 

8.8 

- 2.13 

8.8 

28 

5 

2.92 

8.3 

-.24 

8.2 

28 

6 

- 8.33 

8.1 

2.00 

8.1 

29 

5 

3.45 

7.9 

3.54 

7.9 

29 

6 

-.30 

7.8 

- 2.47 

7.8 

30 

5 

3.39 

7.4 

.70 

7.4 

30 

6 

- 3.22 

7.4 

4.17 

7.4 

31 

5 

.91 

7.1 

1.46 

7.1 

31 

6 

-.56 

7.0 

.77 

7.0 

32 

5 

-.05 

6.7 

- 2.79 

6.7 

32 

6 

- 3.83 

6.6 

.22 

6.6 

33 

5 

-.51 

6.3 

3.24 

6.3 

33 

6 

1.38 

6.4 

- 1.25 

6.3 

34 

5 

- 1.30 

6.0 

.04 

6.0 

34 

6 

.58 

6.0 

-.32 

6.0 

35 

5 

-.24 

5.7 

-.74 

5.7 

35 

6 

.83 

5.7 

- 1.48 

5.7 

36 

5 

- 1.18 

5.4 

.34 

5.4 

36 

6 

-.58 

5.4 

-.90 

5.4 

7 

7 

1.60 

2.1 

22.00 

2.0 







8 

7 

70.42 

1.9 

74.86 

1.9 

8 

8 

- 118.88 

3.9 

122.33 

3.8 

9 

7 

- 118.62 

3.5 

- 100.55 

3.3 

9 

8 

184.50 

3.2 

- 1.85 

3.2 

10 

7 

9.75 

3.4 

- 4.29 

3.5 

10 

8 

43.75 

3.0 

- 92.48 

2.9 

11 

7 

9.61 

4.4 

- 91.89 

4.3 

11 

8 

- 6.35 

4.1 

22.58 

4.0 

12 

7 

- 12.70 

5.2 

34.83 

5.3 

12 

8 

- 21.22 

3.0 

16.90 

3.1 

13 

7 

3.57 

6.5 

- 6.62 

6.3 

13 

8 

- 12.30 

5.0 

- 11.09 

5.0 

14 

7 

37.48 

7.4 

- 4.36 

7.6 

14 

8 

- 32.94 

5.1 

- 13.18 

5.0 

15 

7 

66.71 

8.8 

11.45 

8.8 

15 

8 

- 40.67 

6.7 

24.73 

6.7 

16 

7 

3.05 

10.1 

- 9.07 

10.1 

16 

8 

- 13.44 

8.0 

2.28 

7.9 

17 

7 

22.95 

10.4 

- 11.98 

10.5 

17 

8 

31.16 

9.2 

8.77 

9.2 

18 

7 

-.80 

11.3 

6.72 

11.3 

18 

8 

45.72 

10.8 

.43 

10.5 

19 

7 

5.14 

11.5 

- 1.66 

11.5 

19 

8 

14.86 

10.6 

- 11.33 

10.5 

20 

7 

- 7.79 

11.4 

4.90 

11.4 

20 

8 

- 2.01 

10.8 

- 1.29 

10.8 

21 

7 

- 12.23 

10.8 

- 1.41 

10.7 

21 

8 

- 18.10 

10.7 

2.52 

10.6 

22 

7 

12.75 

11.1 

1.30 

11.1 

22 

8 

- 9.82 

10.2 

- 6.83 

10.1 

23 

7 

- 2.33 

10.4 

2.65 

10.3 

23 

8 

4.22 

9.6 

- 6.80 

9.5 

24 

7 

- 2.51 

10.0 

5.07 

9.9 

24 

8 

- 2.43 

9.7 

7.56 

9.7 

25 

7 

.26 

9.5 

3.46 

9.4 

25 

8 

1.43 

9.1 

- 4.12 

9.1 

26 

7 

5.45 

8.9 

2.58 

8.9 

26 

8 

3.09 

9.1 

- 2.14 

9.1 

27 

7 

6.93 

8.5 

- 2.82 

8.4 

27 

8 

- 4.18 

8.7 

- 4.41 

8.7 

28 

7 

- 4.65 

8.2 

- 1.54 

8.2 

28 

8 

-.57 

8.2 

- 3.18 

8.1 

29 

7 

1.18 

7.9 

- 7.29 

7.8 

29 

8 

- 6.42 

7.7 

2.54 

7.7 




Table 2. Continued. 


1 

m 

Clm 

a 

Sim 

a 

t 

m 

Qim 

a 



a 

~w 

7 

-.17 

7.3 

-.13 

7.3 

30 

8 

2.95 

7.4 

.57 

"TIT 

31 

7 

1.48 

7.0 

- 1.75 

7.0 

31 

8 

.09 

7.0 

- 1.18 

7.0 

32 

7 

- 3.08 

6.7 

1.86 

6.7 

32 

8 

.88 

6.6 

3.32 

6.5 

33 

7 

-.04 

6.3 

1.78 

6.3 

33 

8 

-.08 

6.3 

1.58 

6.3 

34 

7 

2.31 

5.9 

.10 

5.9 

34 

8 

.71 

6.0 

-.79 

6.0 

35 

7 

.07 

5.7 

1.53 

5.7 

35 

8 

.27 

5.7 

.18 

5.7 

36 

7 

-.17 

5,4 

-.42 

5.4 

36 

8 

- 1.05 

5.3 

-.54 

5.3 

9 

9 

- 55.55 

5.9 

97.59 

5.7 







10 

9 

128.18 

3.5 

- 48.19 

3.5 

10 

10 

94.56 

4.3 

- 20.10 

4.2 

11 

9 

- 38.78 

4.1 

40.28 

4.0 

11 

10 

- 52.06 

4.3 

- 17.61 

4.3 

12 

9 

46.94 

4.4 

13.22 

4.5 

12 

10 

- 9.13 

2.8 

31.68 

2.7 

13 

9 

20.38 

4.4 

45.78 

4.4 

13 

10 

43.30 

4.7 

- 38.04 

4.7 

14 

9 

37.16 

5.1 

17.93 

5.2 

14 

10 

37.00 

2.9 

- 2.80 

2.8 

15 

9 

13.44 

5.7 

41.02 

5.7 

15 

10 

9.59 

4.9 

16.08 

4.9 

16 

9 

- 16.58 

6.2 

- 50.98 

6.3 

16 

10 

- 10.42 

4.2 

6.61 

4.1 

17 

9 

- 3.20 

6.9 

- 34.32 

7.0 

17 

10 

2.10 

5.8 

20.12 

5.8 

18 

9 

- 13.52 

8.1 

19.24 

8.1 

18 

10 

9.01 

6.0 

- 10.86 

6.0 

19 

9 

1.76 

8.6 

8.66 

8.6 

19 

10 

- 35.35 

7.4 

- 2.66 

7.3 

20 

9 

22.81 

10.0 

7.24 

10.2 

20 

10 

- 22.42 

8.0 

- 8.09 

8.0 

21 

9 

17.32 

9.8 

- 9.39 

9.8 

21 

10 

3.65 

8.5 

1.84 

8.5 

22 

9 

12.51 

10.2 

- 9.48 

10.3 

22 

10 

5.01 

9.2 

20.38 

9.2 

23 

9 

- 4.03 

10.1 

- 10.39 

10.1 

23 

10 

19.98 

8.8 

- 3.76 

8.8 

24 

9 

- 3.89 

9.8 

- 1.43 

9.8 

24 

10 

17.35 

8.6 

9.30 

8.6 

25 

9 

- 6.03 

9.3 

9.84 

9.3 

25 

10 

5.69 

8.5 

- 4.47 

8.5 

26 

9 

2.52 

8.7 

-.66 

8.8 

26 

10 

- 4.84 

8.5 

1.66 

8.5 

27 

9 

.40 

8.3 

2.16 

8.3 

27 

10 

- 8.31 

8.1 

6.01 

8.1 

28 

9 

2.94 

8.0 

- 3.06 

8.0 

28 

10 

- 7.28 

8.0 

1.20 

8.0 

29 

9 

- 1.64 

7.7 

2.48 

7.7 

29 

10 

.06 

7.5 

6.01 

7.5 

30 

9 

.17 

7.3 

- 3.94 

7.3 

30 

10 

1.23 

7.2 

- 1.04 

7.2 

31 

9 

- 3.89 

6.9 

- 1.70 

6.9 

31 

10 

2.52 

6.8 

- 3.74 

6.8 

32 

9 

1.95 

6.6 

.77 

6.6 

32 

10 

.83 

6.5 

- 1.95 

6.5 

33 

9 

-.29 

6.3 

1.87 

6.3 

33 

10 

.23 

6.1 

-.94 

6.1 

34 

9 

1.27 

5.9 

1.51 

5.9 

34 

10 

- 1.49 

5.9 

-.01 

5.9 

35 

9 

1.08 

5.6 

- 2.07 

5.6 

35 

10 

- 1.43 

5.6 

-.84 

5.6 

36 

9 

-.32 

5.4 

.55 

5.4 

36 

10 

-.36 

5.3 

.47 

5.3 

11 

11 

54.33 

6.2 

- 54.73 

6.3 







12 

11 

5.41 

3.0 

- 9.52 

2.9 

12 

12 

- 3.53 

2.5 

- 11.80 

2.5 

13 

11 

- 40.19 

4.4 

5.50 

4.4 

13 

12 

- 28.01 

2.5 

86.41 

2.5 

14 

11 

8.08 

3.6 

- 41.36 

3.6 

14 

12 

8.97 

1.6 

- 32.07 

1.5 

15 

11 

1.72 

4.3 

28.93 

4.2 

15 

12 

- 28.33 

3.4 

12.49 

3.4 

16 

11 

14.02 

3.9 

- 6.44 

3.9 

16 

12 

20.88 

1.6 

5.74 

1.6 

17 

11 

- 17.11 

4.7 

17.50 

4.7 

17 

12 

34.27 

4.1 

17.26 

4.1 

18 

11 

- 12.80 

4.0 

-.60 

3.9 

18 

12 

- 26.18 

2.4 

- 16.53 

2.3 

19 

11 

16.48 

4.7 

13.47 

4.7 

19 

12 

3.20 

4.7 

4.33 

4.7 

20 

11 

11.38 

4.4 

- 23.93 

4.3 

20 

12 

- 4.06 

3.2 

17.30 

3.1 

21 

11 

9.28 

5.1 

- 36.78 

5.1 

21 

12 

2.82 

5.2 

. 12.71 

5.2 

22 

11 

- 9.37 

5.4 

- 18.38 

5.3 

22 

12 

7.44 

4.0 

- 7.84 

3.9 

23 

11 

3.85 

5.9 

13.68 

6.1 

23 

12 

21.58 

5.5 

- 16.68 

5.6 

24 

11 

12.74 

6.6 

12.12 

6.4 

24 

12 

12.34 

5.0 

- 9.52 

4.9 

25 

11 

5.58 

6.9 

- 1.28 

7.0 

25 

12 

- 5.54 

5.8 

11.01 

5.9 

26 

11 

3.20 

6.8 

5.05 

6.7 

26 

12 

- 19.65 

5.8 

5.43 

5.7 

27 

11 

- 1.19 

7.1 

- 3.10 

7.1 

27 

12 

-.42 

6.1 

- 1.72 

6.2 

28 

11 

-.62 

6.9 

-.83 

7.0 

28 

12 

.40 

6.2 

2.43 

6.2 

29 

11 

- 9.32 

6.9 

.45 

6.8 

29 

12 

-.83 

6.3 

- 4.93 

6.3 

30 

11 

- 1.62 

6.7 

5.08 

6.7 

30 

12 

3.79 

6.2 

- 3.42 

6.2 

31 

11 

- 1.66 

6.4 

5.85 

6.5 

31 

12 

.15 

6.1 

4.62 

6.1 

32 

11 

- 2.40 

6.3 

-.65 

6.3 

32 

12 

- 1.71 

5.8 

4.19 

5,8 

33 

11 

5.56 

6.0 

-.48 

6.1 

33 

12 

5.22 

5.5 

4.13 

5.5 

34 

11 

1.24 

5.7 

- 3.83 

5.7 

34 

12 

.45 

5.3 

2.49 

5.3 

35 

11 

.68 

5.4 

- 3.66 

5.4 

35 

12 

1.48 

5.3 

- 2.06 

5.3 

36 

11 

.59 

5.2 

-.73 

5.3 

36 

12 

-.22 

5.1 

- 1.61 

5.1 

13 

13 

- 61.55 

1.3 

68.27 

1.3 








Table 2. Continued. 


i 

m 

mm 

a 

— Sim 

a 

t 

m 


a 

hh 

a 

14 

13 

31.53 

.9 

44.62 

.9 

14 

14 


1.8 

- 6.37 

1.8 

15 

13 

- 28.11 

1.0 

- 4.98 

1.0 

15 

14 

6.17 

1.0 

- 25.61 

.9 

16 

13 

13.08 

1.1 

.61 

1.1 

16 

14 

- 19.12 

.8 

- 38.29 

.8 

17 

13 

16.91 

1.2 

20.11 

1.2 

17 

14 

- 13.34 

.8 

11.76 

.8 

18 

13 

- 6.58 

1.1 

- 35.16 

1.1 

18 

14 

- 9.28 

1.2 

- 10.94 

1.2 

19 

13 

- 6.09 

1.6 

- 29.17 

1.6 

19 

14 

- 5.12 

.7 

- 12.64 

.7 

20 

13 

26.65 

1.3 

4.89 

1.3 

20 

14 

10.32 

1.8 

- 11.76 

1.8 

21 

13 

- 18.17 

2.1 

11.60 

2.0 

21 

14 

18.78 

1.2 

8.70 

1.2 

22 

13 

- 16.95 

1.4 

17.85 

1.4 

22 

14 

8.73 

2.2 

10.24 

2.2 

23 

13 

- 10.46 

2.8 

- 7.51 

2.8 

23 

14 

4.61 

2.2 

- 3.27 

2.2 

24 

13 

- 3.62 

2.2 

-.38 

2.2 

24 

14 

- 18.64 

2.8 

1.46 

2.8 

25 

13 

7.38 

3.5 

- 15.19 

3.5 

25 

14 

- 21.94 

3.3 

13.21 

3.3 

26 

13 

2.72 

3.2 

1.42 

3.3 

26 

14 

3.93 

3.5 

5.63 

3.5 

27 

13 

- 5.98 

4.2 

- 4.13 

4.1 

27 

14 

11.97 

4.3 

6.64 

4.3 

28 

13 

.10 

3.8 

3.53 

3.9 

28 

14 

- 2.11 

4.3 

- 6.50 

4.3 

29 

13 

- 1.15 

4.5 

- 1.98 

4.5 

29 

14 

- 5.15 

4.6 

1.94 

4.7 

30 

13 

14.67 

4.6 

-.02 

4.6 

30 

14 

-.03 

4.8 

- 2.55 

4.8 

31 

13 

5.69 

4.6 

1.32 

4.6 

31 

14 

- 7.28 

4.4 

1.25 

4.5 

32 

13 

7.26 

5.0 

.22 

5.0 

32 

14 

4.66 

4.6 

6.92 

4.6 

33 

13 

3.67 

4.4 

6.78 

4.3 

33 

14 

9.23 

3.9 

2.51 

3.9 

34 

13 

- 8.08 

4.8 

1.28 

4.9 

34 

14 

- 1.04 

3.9 

-.29 

3.9 

35 

13 

- 1.18 

4.7 

4.48 

4.7 

35 

14 

-.48 

4.2 

-.12 

4.1 

36 

13 

.77 

4.5 

3.79 

4.5 

36 

14 

- 4.84 

4.5 

- 4.07 

4.5 

15 

15 

- 18.09 

3.2 

- 8.09 

3.2 







16 

15 

- 12.53 

3.9 

- 32.30 

3.9 

16 

16 

- 32.41 

7.1 

- 4.37 

7.1 

17 

15 

4.94 

1.4 

5.75 

1.4 

17 

16 

- 29.07 

5.5 

1.88 

5.5 

18 

15 

- 37.76 

3.8 

- 19.82 

3.8 

18 

16 

9.79 

4.0 

5.00 

4.0 

19 

15 

- 18.32 

2.7 

- 12.77 

2.7 

19 

16 

- 19.90 

5.0 

- 11.93 

5.0 

20 

15 

- 22.73 

3.5 

-.41 

3.5 

20 

16 

- 10.67 

4.3 

1.69 

4.3 

21 

15 

16.62 

3.5 

14.98 

3.4 

21 

16 

8.73 

4.2 

- 5.16 

4.2 

22 

15 

27.94 

3.3 

3.10 

3.3 

22 

16 

.09 

5.0 

- 4.93 

5.0 

23 

15 

17.73 

4.1 

- 2.28 

4.0 

23 

16 

4.90 

4.5 

11.77 

4.6 

24 

15 

9.81 

3.6 

- 13.53 

3.6 

24 

16 

-.49 

6.2 

6.28 

6.2 

25 

15 

- 1.99 

4.4 

- 2.27 

4.4 

25 

16 

3.04 

5.3 

- 12.80 

5.3 

26 

15 

- 11.38 

4.8 

4.70 

4.9 

26 

16 

5.82 

6.4 

- 4.15 

6.3 

27 

15 

- 4.34 

5.3 

.10 

5.3 

27 

16 

6.58 

5.5 

- 4.11 

5.6 

28 

15 

- 8.21 

5.8 

5.38 

5.8 

28 

16 

- 8.31 

6.1 

- 7.69 

6.0 

29 

15 

- 1.27 

5.6 

- 2.49 

5.7 

29 

16 

- 2.20 

6.4 

- 5.53 

6.4 

30 

15 

2.81 

6.0 

- 9.28 

6.0 

30 

16 

.62 

6.3 

5.66 

6.3 

31 

15 

.45 

5.6 

- 4.37 

5.6 

31 

16 

- 4.51 

6.4 

4.80 

6.4 

32 

15 

3.91 

5.5 

- 4.94 

5.5 

32 

16 

2.92 

6.2 

4.10 

6.2 

33 

15 

- 3.01 

5.3 

2.17 

5.3 

33 

16 

-.39 

5.7 

1.97 

5.7 

34 

15 

.76 

5.4 

3.00 

5.4 

34 

16 

1.13 

5.6 

- 2.66 

5.6 

35 

15 

.26 

5.1 

2.83 

5.1 

35 

16 

.13 

5.4 

- 1.34 

5.4 

36 

15 

- 1.82 

4.8 

1.86 

4.8 

36 

16 

1.34 

5.2 

- 2.01 

5.2 

17 

17 

- 38.31 

7.5 

- 20.62 

7.6 







18 

17 

6.11 

7.4 

8.77 

7.4 

18 

18 

- 4.45 

11.4 

- 5.06 

11.4 

19 

17 

27.95 

4.0 

- 10.88 

4.0 

19 

18 

21.65 

9.2 

- 3.11 

9.1 

20 

17 

4.29 

7.0 

- 8.98 

7.1 

20 

18 

10.58 

6.5 

1.30 

6.5 

21 

17 

- 6.75 

6.2 

.84 

6.2 

21 

18 

16.83 

8.2 

- 6.57 

8.2 

22 

17 

13.81 

6.4 

- 11.13 

6.4 

22 

18 

7.03 

6.8 

- 10.30 

6.7 

23 

17 

- 7.21 

6.3 

- 6.60 

6.3 

23 

18 

- 1.91 

6.4 

- 6.30 

6.4 

24 

17 

- 8.46 

5.9 

1.88 

5.9 

24 

18 

4.31 

6.5 

- 5.07 

6.5 

25 

17 

- 8.31 

6.4 

.55 

6.4 

25 

18 

- 1.30 

6.7 

- 10.67 

6.6 

26 

17 

- 4.89 

5.9 

8.28 

5.9 

26 

18 

- 9.01 

7.2 

7.55 

7.3 

27 

17 

5.55 

6.4 

1.59 

6.4 

27 

18 

- 5.19 

6.9 

5.94 

6.9 

28 

17 

4.52 

6.9 

- 4.26 

6.9 

28 

18 

.36 

6.7 

-.83 

6.6 

29 

17 

4.57 

6.5 

- 2.78 

6.5 

29 

18 

- 2.01 

6.3 

-.14 

6.3 

30 

17 

1.01 

6.5 

- 1.54 

6.5 

30 

18 

.30 

6.5 

.85 

6.5 

31 

17 

- 5.93 

6.4 

2.57 

6.4 

31 

18 

2.56 

6.4 

.85 

6.4 

32 

17 

- 3.66 

6.3 

1.89 

6.3 

32 

18 

2.21 

6.2 

- 1.48 

6.2 

33 

17 

- 2.14 

6.0 

3.04 

6.0 

33 

18 

.87 

6.1 

- 1.24 

6.1 



Table 2. Continued. 
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17 

.36 
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2.58 

5.7 

34 

18 

- 2.72 
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-.09 

5.9 

35 

17 

3.39 

5.5 

- 2.44 

5.5 

35 

18 

1.50 

5.5 

-.59 

5.5 

36 

17 

2.17 

5.3 

-.85 

5.3 

36 

18 

.07 

5.2 

.53 

5.2 

19 

19 

6.46 

10.5 

10.42 
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20 

19 

- 7.10 

9.4 

8.46 

9.5 

20 

20 

1.71 

13.1 

- 13.51 

13.2 

21 

19 

- 20.95 

5.7 

15.88 

5.7 

21 

20 

- 19.04 

9.9 

18.54 

9.9 

22 

19 

6.62 

8.2 

- 4.70 

8.3 

22 

20 

- 13.32 

7.6 

14.78 

7.6 

23 

19 

- 8.68 

7.9 

7.49 

8.1 

23 

20 

17.22 

8.4 

- 9.05 

8.5 

24 

19 

.53 

7.3 

- 15.02 

7.2 

24 

20 

- 6.06 

8.0 

-.33 

8.1 

25 

19 

9.18 

6.9 

2.13 

7.0 

25 

20 

- 3.73 

7.0 

- 6.62 

6.9 

26 

19 

1.63 

6.7 

.72 

6.7 

26 

20 

9.48 

6.8 

- 10.95 

6.8 

27 

19 

.93 

6.7 

- 6.23 

6.7 

27 

20 

2.95 

6.7 

3.01 

6.7 

28 

19 

4.42 

6.7 

13.80 

6.7 

28 

20 

-.98 

7.0 

1.13 

7.0 

29 

19 

- 2.19 

6.8 

1.51 

6.8 

29 

20 

- 4.80 

6.7 

3.16 

6.7 

30 

19 

- 5.63 

6.7 

- 2.64 

6.7 

30 

20 

-.03 

6.5 

3.57 

6.5 

31 

19 

2.02 

6.2 

2.82 

6.3 

31 

20 

1.89 

6.3 

.56 

6.2 

32 

19 

3.09 

6.2 

- 1.51 

6.2 

32 

20 

- 1.59 

6.2 

1.42 

6.2 

33 

19 

1.68 

5.9 

.24 

6.0 

33 

20 

2.08 

5.9 

-.75 

5.9 

34 

19 

.60 

5.7 

- 1.00 

5.7 

34 

20 

.89 

5.8 

-.40 

5.7 

35 

19 
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5.6 

.59 

5.6 

35 

20 

-.77 

5.5 

-.88 

5.6 

36 

19 

-.26 

5.3 

-.07 

5.3 

36 

20 

-.91 

5.3 

-.87 

5.3 

21 

21 

2.48 

11.3 

- 6.85 

11.3 
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21 

— 13.22 

10.4 

7.60 

10.4 

22 

22 

- 1.46 

12.7 

4.72 

12.6 

23 

21 

10.82 

6.4 

7.64 

6.4 

23 

22 

-.90 

9.5 

- 2.14 

9.5 

24 

21 

10.57 

8.2 

1.12 

8.1 

24 

22 

- 1.73 

7.9 

- 1.30 

7.9 

25 

21 

5.40 

8.0 

3.14 

8.0 

25 

22 

- 1.87 

7.5 

- 1.76 

7.5 

26 

21 

-.39 

7.3 

- 2.41 

7.3 

26 

22 

10.91 

7.8 

9.16 

7.8 

27 

21 

2.07 

6.7 

- 4.54 

6.6 

27 

22 

-.14 

6.6 

2.93 

6.6 

28 

21 

2.49 

6.6 

.27 

6.6 

28 

22 

- 4.85 

6.5 

.51 

6.4 

29 

21 

- 9.40 

6.2 

- 5-95 

6.1 

29 

22 

9.66 

6.3 

4.44 

6.3 

30 

21 

- 7.45 

6.5 

- 3.10 

6.5 

30 

22 

3.25 

6.3 

- 5.54 

6.3 

31 

21 

2.34 

6.4 

3.58 

6.4 

31 

22 

- 6.27 

6.0 

- 5.73 

6.0 

32 

21 

1.17 

6.1 

5.88 

6.1 

32 

22 

- 4.66 

6.1 

.75 

6.1 

33 

21 

.73 

5.9 

-.84 

5.8 

33 

22 

- 4.10 

5.8 

- 1.10 

5.8 

34 

21 

1.37 

5.7 

-.66 

5.8 

34 

22 

.83 

5.6 

.42 

5.6 

35 

21 

1.38 

5.4 

2.74 

5.4 

35 

22 

.02 

5.4 

3.30 

5.3 

36 

21 

.74 

5.2 

- 2.12 

5.2 

36 

22 

.57 

5.2 

.69 

5.2 

23 

23 

.84 

10.5 

.20 

10.5 







24 

23 

- 2.14 

9.9 

- 9.01 

9.9 

24 

24 

2.34 

11.2 

- 1.21 

11.2 

25 

23 

4.57 

6.3 

- 2.46 

6.3 

25 

24 

3.61 

8.5 

- 3.86 

8.5 

26 

23 

2.35 

7.1 

8.95 

7.1 

26 

24 

- 1.37 

6.9 

12.18 

7.0 

27 

23 

- 5.39 

7.0 

- 2.74 

7.0 

27 

24 

- 1.94 

6.1 

2.62 

6.1 

28 

23 

- 2.64 

6.7 

6.39 

6.6 

28 

24 

6.88 

6.3 

- 15.08 

6.3 

29 

23 

- 5.01 

6.0 

-.08 

6.0 

29 

24 

- 2.53 

5.4 

3.70 

5.4 

30 

23 

- 1.58 

6.0 

- 5.34 

5.9 

30 

24 

- 2.51 

5.6 

-.04 

5.6 

31 

23 

9.50 

5.4 

5.65 

5.5 

31 

24 

- 3.83 

5.5 

- 1.91 

5.4 

32 

23 

3.86 

5.8 

.48 

5.7 

32 

24 

- 6.59 

5.3 

5.33 

5.3 

33 

23 

-.63 

5.6 

- 4.40 

5.6 

33 

24 

3.99 

4.9 

-.48 

4.9 

34 

23 

-.93 

5.5 

- 2.14 

5.5 

34 

24 

6.73 

5.3 

.85 

5.3 

35 

23 

- 2.38 

5.2 

- 1.56 

5.2 

35 

24 

2.53 

5.0 

2.20 

5.0 

36 

23 

- 1.22 

5.0 

-.55 

5.0 

36 

24 

.66 

4.9 

- 1.42 

4.9 

25 

25 

4.95 

8.5 

4.01 

8.5 







26 

25 

- 3.99 

8.3 

8.25 

8.4 

26 

26 

3.43 

7.6 

- 4.27 

7.5 

27 

25 

11.80 

5.5 

3.15 

5.5 

27 

26 

- 5.01 

6.4 

4.00 

6.4 

28 

25 

1.13 

6.4 

- 4.81 

6.5 

28 

26 

3.44 

5.7 

1.68 

5.7 

29 

25 

8.31 

6.1 

3.63 

6.0 

29 

26 

6.26 

5.0 

- 3.69 

5.1 

30 

25 

8.85 

6.0 

- 5.61 

6.0 

30 

26 

- 3.21 

5.9 

8.14 

5.9 

31 

25 

- 7.74 

5.2 

-.28 

5.2 

31 

26 

- 4.69 

4.8 

.36 

4.8 

32 

25 

- 13.17 

5.1 

7.70 

5.1 

32 

26 

- 1.05 

4.9 

- 1.53 

4.9 

33 

25 

- 1.26 

4.7 

- 4.82 

4.7 

33 

26 

8.13 

4.3 

5.54 

4.3 

34 

25 

6.25 

4.9 

- 8.21 

4.9 

34 

26 

1.08 

4.6 

- 9.01 

4.6 

35 

25 

- 3.39 

4.7 

1.52 

4.7 

35 

26 

- 14.34 

3.7 

-.14 

3.7 
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Qssl 

ie 25 j03 

27 27 6.89 

27 - 9.92 

27 - 7.47 

27 - 1.92 

27 7.00 

27 - 3.09 

27 - 10.33 

27 6.90 

27 2.78 

27 - 10.19 

29 8.64 

29 4.82 

29 - 5.42 

29 - 3.34 

29 - 21.32 

29 - 3.81 

29 8.19 

29 - 1.34 

31 -.23 

31 -.76 

31 .17 

31 2.35 

31 1.33 

31 - 3.55 
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33 - 2.40 
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4.9 
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1.33 
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4.7 
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4.7 

3.5 
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4.4 
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4.4 

2.6 
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2.6 

4.0 

-.54 
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4.23 

2.6 
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3.16 
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6.3 
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3.9 
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4.0 

4.1 

2.50 

4.0 

3.5 

.19 

3.5 

3.1 

- 4.48 
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3.3 

- 3.15 

3.3 

2.9 

- 2.45 

2.9 

7.2 

-.07 

7.2 

6.3 

- 2.71 

6.3 

5.6 

1.16 

5.6 

4.3 

2.30 

4.3 

5.1 

1.03 

5.1 

4.9 

1.36 

4.8 

6.5 

-.36 

6.5 

6.0 

1.38 

6.0 

5.6 

1.67 

5.6 

4.9 

- 2.62 

4.9 

5.8 

-.14 

5.8 

5.4 

-.73 

5.4 
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26 
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28 

28 

6.77 

29 

28 

10.32 

30 

28 
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31 
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.34 

32 

28 

1.57 

33 

28 
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34 

28 

4.55 

35 

28 
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36 

28 

6.95 

30 

30 
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31 

30 

- 2.45 

32 

30 

8.29 

33 

30 
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34 

30 
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35 

30 

3.71 
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32 

32 

.73 

33 

32 
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35 

32 
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3.9 

4.7 
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4.7 
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4.1 
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3.3 
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4.4 
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6.8 

6.3 
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6.3 

5.4 

- 1.16 

5.4 

5.3 
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6.1 

-.60 

6.1 

5.7 

.12 

5.7 

5.2 

1.76 

5.2 


36 36 
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3 Satellite Dynamics 


3.1 Cartesian Equations of Motion 

The gravitational acceleration at any given location is obtained by computing the 
gradient of the potential. Since the potential is given as a function of Earth-fixed 
spherical coordinates, it is most convenient to compute the gradient in the same 
system. In Earth-fixed spherical coordinates, this gradient is, 


a = VC/ = 


dU ^ 1 dU ~ 1 dU ~ 

—Ur + -r-u^ H T 

dr r d<t> r cos 4> d\ 


(13) 


where u r , and Hx are unit vectors in the r, 4>, A basis. This basis has u P pointing 
along the radius vector to the satellite, u^ is in the direction of increasing north lati- 
tude and u\ is in the direction of increasing east longitude. The acceleration vector 
obtained from this expression will be the inertial acceleration for the point of interest. 
Though, as noted, the components of the acceleration are given in the Earth-fixed 
coordinate system. For most applications it will be desired to have the components of 
the acceleration expressed in an inertial (nonrotating) coordinate system. This is ac- 
complished by applying the appropriate coordinate transformation from the spherical 
coordinates to the desired coordinate system. So as a first step, the components of 
the inertial acceleration in the Earth-fixed (rotating) coordinate system are obtained. 
Substituting for the gravitational potential (equation (1)) and taking the indicated 
partials in equation (13) gives the acceleration vector, 


a = 1) ) E P em (sm<t>)[C tm cosm\ + St m sinm\]\u r 

l r e=o \ r / m _o J 


{55 (t)'s 

{S£(?)'S 


M* XT' dPfcn (sin <(>) 


d<t> 


[Ce m cos mA -I- S( m sin mA] > 


1 


m fim(sin 0) y_c tm si n m \ + S (m cos mA] ^ u x (14) 
cos <fi 1 


I) 


Notice that the leading term of the radial component (degree and order equal to zero) 
is simply the expected two-body gravitational acceleration —n/r 2 . Also, if only zonal 
terms are used (m = 0), then the longitudinal component of the acceleration is zero. 

Next, the Earth-fixed Cartesian components of the acceleration can be obtained 
by rotating from the spherical coordinates to the x, y, z basis. Let the components 
of the acceleration in spherical coordinates be represented by, 


a — a r u r + a 0 u* + a x ux 


( 15 ) 


where the components a r , a# and a* are given in equation (14). The acceleration 
vector in Cartesian coordinates can be written as, 


Clxyz — UjUj "f" UyUy “I" fljUj 


( 16 ) 



where u x , v„ and u z are the Cartesian unit vectors in the Earth-fixed (rotating) 
coordinate system. The Cartesian components of the acceleration can be obtained 
from the spherical coordinate components through the standard transformation, 

a x cos 0 cos A — sin0cosA —sin A a r 

a v = cos0sinA — sin0sinA cos A (17) 

a, sin0 cos0 0 a\ 

Having obtained the Earth-fixed Cartesian components of the acceleration one 
further coordinate transformation is necessary to obtain the acceleration components 
in the defined inertial coordinate system. If the matrix T represents the coordinate 
transformation from the Earth-fixed system to the inertial coordinate system, then 
the acceleration components in the inertial system will be, 

Qxyz — Ta X yz (IS) 

where a XYZ is the inertial acceleration vector in inertial coordinates, 

&XYZ — Qx u x flyUy + Q Z tl Z (19) 

with u x , tTy and u z being the unit vectors of the Cartesian inertial coordinate system. 
In component form, this final transformation will have the structure, 


ax 

Tn 

Tn 

Tn a* 


Qy 

= Tn 

Tn 

T23 a y 

(20) 

a z 

Tn 

Tyz 

T33 a z 



The actual elements of the transformation matrix T depend on the inertial coor- 
dinate system being used. In the most general case, this transformation will account 
for polar motion (the motion of the spin axis with respect to the Earth crust), Earth 
rotation (the largest effect) and, precession and nutation (the motion of the spin axis 
with respect to the stars). In the simplest case, all of these effects are neglected 
except for Earth rotation. This defines a coordinate system with the same z axis as 
the Earth-fixed system but not rotating with the Earth. For many applications such 
a system is effectively inertial. The transformation from the Earth-fixed system to 
this nonrotating system is simply, 

cos 8 — sin 8 0 

T = sin# cos 8 0 (21) 

0 0 1 

where 8 is the Greenwich Hour Angle (the angle from a reference direction, usually 
the Vernal Equinox, to the Greenwich meridian). 



3.3 Lagrangian Equations of Motion 


The equations of motion for a satellite moving in the gravity field are given by equation 
(13). These equations are convenient for numerical computation of an ephemeris but 
not for analytic investigations into the evolution of the orbital elements. Such analyses 
are more conveniently performed using the Lagrangian equations of motion. These 
equations are equivalent to equation (13) but directly give the time rates of change 
of the Kepler orbital elements, 


da 


2 dV 

(22) 

dt 

— 

na dM 

de 

dt 
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na 2 e dM na, 2 e du 

(23) 
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(24) 
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(25) 

dt 
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dQ 


1 dV 

(26) 

dt 


na 2 (l — e 2 ) 1 / 2 sini di 

dM 


1 -e 2 dV 2 dV 

(27) 

dt 

= 

na 2 e de na da 




(28) 


where V is the perturbing gravitational potential, i.e., the gravitational potential 
excluding the point mass contribution, 


U = - + V (29) 

r 

These equations require that the perturbing gravitational potential V be repre- 
sented in terms of the orbital elements instead of the spherical coordinates as ex- 
pressed in equation (1). This conversion has been carried out by Kaula [1966], 


“ /=! ' “ 7 msOisO <?=— oo 


where 


and 


Stmpq — 


c (m 

~S( m 


I t—m even 


COS Iptmpq + 


I t—m odd 


S( m 


t—m even 


sin Vy, 


mpq 


J t—m odd 


4>tmpq = (/ - 2p)u + (/ - 2p + q)M + m(Cl - 9) 


(30) 

(31) 

(32) 


All quantities having been previously defined except for the inclination function 
Ftmp(i) and eccentricity function Gt„{e), which are given by Kaula. In this for- 
mulation, the overbar on F( mp indicates it is normalized (in the same fashion that 



3.2 Orbital Elements 


While it is convenient to express the satellite equations of motion in the Cartesian 
coordinate system, the evolution of the satellite orbit is more conveniently evaluated 
in terms of its orbital elements. This is because the satellite motion is nearly Ke- 
plerian (i.e., it corresponds closely to the ideal elliptical motion which results from 
the gravitational motion about a point mass). The orbital elements are a set of six 
parameters which uniquely define the position and velocity of the satellite. Transfor- 
mations exist which convert Cartesian position and velocity to orbital elements and 
vice-versa. The set of elements used in this study are the conventional set of Kepler 
elements. These six elements are 

semi-major axis — a — One-half the length of the major axis of the ellipse. For a 
circular orbit this would be the radius of the circle. 

eccentricity — e — Measure of ellipticity (0 < e < 1). For a circular orbit e = 0. 

inclination — i — The angle formed by the orbital plane and the equatorial (x-y) 
plane (0° < i < 180°). 

right ascension of the ascending node — Q — The angle measured in the equa- 
torial (x-y) plane from the x-axis to the point of intersection of the orbital path 
with the equatorial plane. The orbital path will intersect the x-y plane in two 
places: once ascending (south to north) and once descending (north to south). 
The x-axis will normally be defined to be pointing at the vernal equinox (the 
location of the Sun in the sky on or about March 21). 

argument of perigee — u — The angle measured in the plane of the orbit from 
the ascending node to the point of the orbit corresponding to closest approach 
(the distance from the center of the Earth to the satellite is a minimum). This 
angle will not be defined if the orbit is circular. 

mean anomaly — M — A non-physical angle. Defined as A/ = n(t — T) where T 
is the time of periapsis passage, t is the current time, and n is the mean motion 
(n = 2 tt/P, where P is the orbital period). Unless the orbit is circular, the 
satellite does not move with a constant angular rate and thus the actual angle 
swept out by the satellite (which is the true anomaly) in a given amount of 
time will not correspond to the mean anomaly. The mean anomaly would be 
the angle swept out by the satellite if it were moving at a uniform rate. 

At any given instant the Cartesian position and velocity of the satellite can be con- 
verted to the corresponding set of Kepler orbital elements. These values of the or- 
bital elements will continually change in time due to the various forces affecting the 
satellite; these instantaneous values for the orbital elements are referred to as the 
osculating orbital elements. And it is the changes in these osculating orbital elements 
(due to the acting forces) that is of primary interest here. 



the Associated Legendre Functions are normalized via equation (9)). This is required 
since normalized spherical harmonic coefficients are being utilized. 




4 Satellite Motion 


4.1 Precessing Ellipse 


A satellite moving in the low-Earth orbit environment does not follow the idealized 
elliptical path resulting from the consideration of the motion of a particle about a 
point mass. Instead, the satellite will deviate from the idealized elliptical path as 
the various acting forces perturb the satellite. In this discussion, only the perturbing 
forces resulting from the Earth’s gravity field are considered. Though it should be 
emphasized that depending on the actual orbital geometry and spacecraft design, 
other forces (principally, atmospheric drag, solar radiation pressure, and lunar and 
solar gravity) may be more important. 

The general problem of determining the motion of a satellite orbiting in a gravita- 
tional field (described by spherical harmonics) can only be solved through numerical 
integration. The equations presented in section 3 can be used for this purpose. How- 
ever, it is possible through various degrees of approximation to analytically determine 
the important characteristics of the motion. 

The largest effect on a satellite orbiting the Earth is that due to the oblateness 
(predominantly the second degree zonal C^o). This is not surprising given the size 
of the second degree zonal with respect to all other coefficients of the gravity field 
(Tables 1 and 2). The most noticeable effect of the oblateness is to cause the right 
ascension of the ascending node and argument of periapse of the orbit to vary linearly 
in time. In the two-body problem these quantities remain constant. The oblateness 
also causes a small change in the orbital period. Besides these secular effects (linear 
changes in time) the oblateness causes (relatively) large periodic changes in all of the 
orbital elements. 

The rate of change in the node, periapse, and mean anomaly due to the second 
degree zonal coefficient are (from Kaula [1966]), respectively, 


dfi 3 /a e \ 2 cosi 

It ~ 2 n Uj (1 - e 2)2 C2 - 0 

du _ 3_ ( a e \ 2 (5cos 2 i — 1)^, 

Tt = TUJ c ™ 

dM 3 /a e \ 2 (3cos 2 i — 1) 

IT = n "4 n uJ (l-e 2 )^ C2 ’° 


(33) 

(34) 

(35) 


(Note that the unnormalized second degree zonal coefficient is used in these expres- 
sions.) 

So as a first approximation, a satellite orbiting the Earth can be reasonably rep- 
resented by a secularly precessing ellipse. This is an ellipse that has constant values 
for the semi-major axis, eccentricity and inclination, and values that vary linearly in 
time for the right ascension of the ascending node, argument of periapse and mean 
anomaly. The actual motion of the satellite then consists of small deviations away 
from this secularly precessing ellipse. To characterize these deviations away from this 
precessing ellipse the solution presented by Kaula [1966] is very useful. 



4.2 Kaula’s Solution 

Kaula’s approach was to define a reference orbital path corresponding to the secularly 
precessing ellipse. The orbital elements that define this reference orbital path are 
referred to as mean (or averaged) orbital elements. Kaula determined (to first order) 
the deviation from these mean orbital elements caused by the spherical harmonic 
coefficients of the gravity field. That is, it is assumed that the orbital elements can 
be written in the form, 


semi-major axis 

a* = 

a + A a 

(36) 

eccentricity 

e m = 

e + Ae 

(37) 

inclination 

i* = 

i + Ai 

(38) 

right ascension of the ascending node 

ft- = 

ftt + ft 0 4" Aft 

(39) 

argument of periapse 

w* = 

tit + u„ + A u 

(40) 

mean anomaly 

M* = 

Mt + M 0 + AM 

(41) 


where the star quantities represent the osculating (true) orbital elements of the satel- 
lite and the unstarred orbital elements are those that define the reference precessing 
ellipse. The time t is measured from a reference epoch and ft<„ and M 0 are the 
values at the reference epoch. The deviations from the reference ellipse will be purely 
periodic (assuming the mean orbital elements and rates have been accurately defined). 
The resulting deviations in the orbital elements are of the form 

Ao = f:EEE Aa, m „ (42) 

{sal m*0 p=0 <7=-oc 

where a represents any one of the six orbital elements. The individual Aa< mw can 
be written in the forms 

A a lm „ = Cf mpq S tm „ for a, e, and i (43) 

A a tm „ = C? mp9 S m tmpq for ft, w, and M (44) 

where the Cf mpq are constants which depend upon the orbital elements of the reference 
secularly precessing ellipse; and the St mpq and S( mpq functions are sinusoidal with 
amplitudes dependent on the values of the gravity field spherical harmonic coefficients. 
The constant Cf mpq coefficients are 

r '« 
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C* 
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The periodicities of the orbital element variations caused by the gravity field are 
thus, 

i>tm„ = (£- 2 p)u + (£-2p + q)M + m(fi - 9) (53) 

where 6 is the Earth rotation rate, £ and m are the spherical harmonic degree and 
order, p is an integer that can take on the values 0 to £, and q is an integer that can 
take on the values ±oo. Kaula’s solution is obtained through an expansion in the 
orbital eccentricity which is the reason the index q has an infinite extent. For orbits 
with small eccentricities only a small range of q values needs to be considered. High 
eccentricity orbits require the consideration of much greater range in q and Kaula’s 
solution may not be the best approach in such cases. 


4.3 Position Perturbation 

Kaula’s solution provides a prediction of the expected variation in the Kepler or- 
bital elements due to the gravity field. This formulation is useful for determining 
the amplitudes and frequencies of the variations in the osculating orbital elements. 
To determine how these orbital element perturbations are manifested in terms of a 
position displacement of the satellite (with respect to the reference secularly process- 
ing ellipse) it is necessary to transform the results back into a Cartesian coordinate 
system. The chosen coordinate system for doing this is the rotating radial, transverse 
(along-track), and normal (cross-track) system. This coordinate system is defined by 
the reference secularly precessing ellipse. The radial direction is along the radius vec- 
tor, tr ans verse is perpendicular to the radius vector (in the orbital plane and positive 
in the direction of the satellite motion), and normal is perpendicular to the orbital 
plane (and positive in the direction of the angular momentum vector). 


The corresponding position perturbation in the radial, transverse, and normal 
directions is given by the following relations (Rosborough and Tapley, [1987]) 


Ar = 

dr dr dr 

-Aa + + —AM 

(54) 

Ar = 

r (au> + + Aftcosi) 

(55) 

A 7 ] = 

r [At sin(u/ + /) — A Cl sin t cos(u> + /)] 

(56) 


where / is the true anomaly as measured along the reference secularly processing 
ellipse (the central angle from the point of periapse to the satellite location). Given 
the expected perturbations in the orbital elements from Kaula’s solution it is then 
possible to compute the corresponding position perturbations using these mapping 
relations. 

In the application at hand, the satellites of interest are established in orbits with 
small eccentricity. This allows for an accurate evaluation of the expected orbital el- 
ement variations using a very limited range of q values. This in turn simplifies the 
evaluation for the position perturbation. Assuming a small eccentricity, the pertur- 
bations in the radial, transverse and normal components are approximately given 
by 


Ar = 21 5Z H CimpStmp 0 ( 57 ) 

1= 1 m=0p*0 

Ar = EIEVw < 58 > 

U1 m=0 p=»0 

Ar/ = (ptmp^{t+l)mpO — ^tmp^(t-l)mpo) (59) 

£=1 m=0p=0 

The constant coefficients C£ mp , C7 mp , C^ p and Cf~ p depend upon the orbital elements 
and rates of the reference secularly processing ellipse and are given by, 
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5 Perturbations in Laser Range 

Of particular interest in this application is the expected orbit error signal that will 
be observed using satellite laser ranging. The orbit perturbations produced by the 
gravity field cause a perturbation in position (in three-dimensions) as described in the 
preceding section. This perturbation in position will result in a perturbation of the 
expected station- to-satellite range measurement. The final step is to then determine 
this perturbation in the range observable. 

The station- to-satellite range is 

P = \J[x- x,) 2 + (y - y,) 2 + (z - z,) 2 (64) 

where x, y, z are the Cartesian coordinates of the satellite and x„ y,, z, are the Carte- 
sian coordinates of the tracking station. The perturbation in the range measurement 
due to a perturbation in the satellite position is given (to first order) by 


A p 


a? p — v-y« £r.£x 

A a p p p 


Ax 

Ay 
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C J. Ax + lZ«± A y + i-±Z Z 
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(65) 

( 66 ) 


where Ax, Ay, and A z axe the perturbations in the satellite position. The perturba- 
tions in the satellite Cartesian position will now be expressed in terms of the already 
specified radial, transverse, and norm ad position perturbations. The transformation 
between these two orthogonal coordinate systems is 
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The rows of the R rotation matrix contain the components of the unit vectors in the 
radial, transverse, and normad directions. That is, 
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U r 
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U v x Ur 
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( 68 ) 


Thus, given the perturbations in the radial, tramsverse, and normal components, it 
is possible to compute the corresponding perturbations in the xyz Cartesian system. 
And the resulting perturbation in the range measurement is 


A p = A t B = A t R t D = A t R t 


Ar 

At 

At? 


( 69 ) 



6 Covariance Analysis 

The preceding analysis provides the tools for determining the perturbations in the 
satellite orbit position (in terms of orbital elements or Cartesian position) and the 
resultant perturbation in a station-to-satellite range measurement. These perturba- 
tions describe the effect on the orbit of the perturbing gravitational potential with 
respect to the reference precessing ellipse. Since the relations are all linear with re- 
spect to the gravity spherical harmonic coefficients, one can equally well determine 
the effect of errors in the coefficients (as opposed to the entire effect). Since a nominal 
gravity model is used in the analysis of geodetic satellite tracking data, it is in fact 
the error in the coefficients that is of most interest and not the entire effect (most of 
which is being accurately modeled by the nominal gravity model). However, the error 
in the coefficients is obviously not known. But the statistics of the expected error 
(in the gravity field coefficients) is known. This knowledge of the error statistics is a 
by-product of the least-squares process used in the estimation of the nominal gravity 
model, and it is represented by the estimated error covariance matrix. 

Define a vector which contains the differences of the true gravity field spherical 
harmonic coefficients with respect to the nominal gravity model values 

8C20 

8C 21 

5521 

8O& 

5522 

8Cm (70) 

8C 3 , 

55 3 i 

5C40 


where the coefficients with an asterisk denote the values used in the nominal gravity 
field model. The error covariance matrix of the estimated coefficients is then 

P g = E[8gSg T ] (71) 

where E is the expectation operator. The gravity field coefficient error covariance 
matrix P g is an n x n symmetric matrix where n is the total number of coefficients. 

Using the gravity field error covariance and the linear relations between the orbit 
perturbations and the gravity field coefficients, it is possible to compute the estimated 
orbit position variance. Similarly, the estimated station-to-satellite range error vari- 
ance can also be obtained. So for a given uncertainty in the gravity field (as provided 
by the estimated error covariance matrix P g ) it is possible to determine what the 
expected uncertainty in the satellite position and/or range measurement will be. 


C20 — Qo 
C21-Q1 

521 — 521 
C22 — 0^2 

522 ~ Sjl2 

8g= 5 30 “ 

On — On 

5n ~ $n 

C40 — C40 


As noted, the relations between the orbit perturbations and the gravity field co- 
efficients are linear, and thus can all be expressed in the form 

6e = F T 6g (72) 

where F is the vector of coefficients that multiply the gravity field coefficients. The 
quantity Se can represent any quantity of interest that has been previously derived: 
orbital element; radial, transverse, normal position component; Cartesian position 
components; or the range measurement error. In each case, the coefficients in the 
vector F will differ, but the linear relation with respect to the gravity field coefficients 
always exist. The lower case <5 is used to indicate that the perturbation is resulting 
from an error in the gravity field coefficients; as opposed to the previously used upper 
case A which was used to denote the full perturbation effect due to the gravity field 
coefficients. 

The variance in the quantity e is then 

P ( = E[F T 5gSg T F] = F T E[SgSg T )F = F T P g F 


(73) 


7 Applications 

These analytical developments were applied to the particular case of laser range track- 
ing of the Lageos satellite. The problem was to determine the expected magnitude 
and characteristic of the Lageos orbit error due to gravitational field mismodeling 
as observed through laser range tracking. In this study the gravity field error co- 
variance for GEM-Tl and GEM-T2 were used to quantify the current uncertainty 
of the Earth’s gravity field. Thus the results presented here are only indicative of 
the current state of the art as represented by GEM-Tl and GEM-T2. Though the 
methodology described in this report can continually be applied to new fields that 
are developed and thus provide another means for quantifying the improved accuracy 
that such fields provide. 

7.1 Lageos and the Tracking Network 

The geodetic satellite Lageos is am excellent satellite for analysing the effect of gravity 
field mismodeling. The satellite is in a high near-circular orbit which decreases its 
sensitivity to the high degree gravity field coefficients but at the same time it is 
also minimally effected by nongravitational forces which could mask the gravitational 
perturbations. So while it does not provide a good test for the entire gravity field 
model, it does provide a stringent test for the low degree portions of the gravity field 
(approximately up to degree 20). The orbitail elements of the Lageos satellite are 
given in Table 3. 


Table 3. Lageos Orbitad Elements. 


Orbital Element 

Value 

semi-major axis 

12,271 km 

eccentricity 

0.0044 

inclination 

109.84 degrees 


Geodetic satellites such as Lageos are tracked by a globally distributed network of 
laser ramging stations. Some of these stations operate at a fixed location and others 
are tramsportable amd have operated at a number of locations. This has resulted in 
observations being obtained at well more than 50 locations. To evaluate the expected 
Lageos laser range residuals that could result from mismodeling the gravity field, a 
representative subset of 20 tracking locations have been used in this study. These 
stations are globally distributed and are listed in Table 4. 



Table 4. Subset of Laser Range Tracking Station Locations. 


Location 

East Longitude 

Latitude 

Arequipa, Peru 

288.5° 

-16.5° 

Easter Island, Chile 

250.6° 

-27.1° 

Goldstone, USA 

243.2° 

35.2° 

Grasse, France 

6.9° 

43.8° 

Graz, Austria 

15.5° 

47.1° 

Greenbelt, USA 

283.2° 

39.0° 

Haleakala, USA 

203.7° 

20.7° 

Huahine, French Polynesia 

209.0° 

-16.7° 

Kwajalein, USA 

192.5° 

9.4° 

Matera, Italy 

16.7° 

40.6° 

McDonald Observatory, USA 

256.0° 

30.7° 

Monument Peak, USA 

243.6° 

32.9° 

Orroral, Australia 

149.0° 

-35.6° 

Platteville, USA 

255.3° 

40.2° 

Quincy, USA 

239.1° 

40.0° 

Royal Greenwich Observatory, UK 

0.3° 

50.9° 

Shanghai, China 

121.2° 

31.1° 

Simosato, Japan 

135.9° 

33.6° 

Wettzell, Germany 

12.9° 

49.1° 

Yaragadee, Australia 

115.3° 

-29.0° 


7.2 Predicted Range Uncertainty 

For each of the tracking station locations listed in Table 4, the expected uncertainty 
of the station-to-satellite range was determined for the Lageos satellite and using the 
error covariance matrices of the GEM-T1 and GEM-T2 gravity field models. These 
uncertainties were computed for every geographic location of the Lageos satellite that 
placed the satellite at least 20 degrees above the horizon with respect to the tracking 
station location. Thus, the uncertainty varied depending on the relative orientation 
of the satellite with respect to the tracking site. The results of these computations 
are summarized in Table 5 for GEM-T1 and Table 6 for GEM-T2. 


Table 5. Station-to-Satellite Range Uncertainty for Lageos. 


Based on the GEM-T1 Gravity Field Error Covariance. 




Range Uncertainty 

in millimeters 



Ascending 


Descending 

Overall 

Station 

Min 

Max RMS 

Min 

Max 

RMS 

RMS 

Arequipa 

11 

22 

16 

13 

23 

17 

16 

Easter Island 

12 

33 

20 

12 

24 

17 

19 

Goldstone 

11 

25 

17 

12 

20 

15 

16 

Grasse 

11 

26 

18 

12 

33 

22 

20 

Graz 

11 

28 

18 

12 

33 

22 

20 

Greenbelt 

11 

23 

16 

12 

23 

16 

16 

Haleakala 

14 

31 

21 

12 

21 

16 

18 

Huahine 

15 

35 

23 

12 

23 

17 

20 

Kwajalein 

15 

33 

22 

12 

21 

16 

19 

Matera 

12 

24 

18 

12 

35 

22 

20 

McDonald Observatory 

11 

23 

16 

12 

21 

15 

16 

Monument Peak 

11 

24 

17 

12 

20 

15 

16 

Orroral 

12 

32 

21 

11 

24 

17 

19 

Platteville 

11 

23 

17 

11 

21 

16 

16 

Quincy 

11 

25 

18 

11 

21 

15 

17 

Greenwich Observatory 

11 

27 

18 

11 

30 

21 

19 

Shanghai 

13 

31 

21 

13 

28 

19 

20 

Simosato 

13 

33 

22 

12 

26 

18 

20 

Wettzell 

11 

29 

18 

12 

33 

22 

20 

Yaragadee 

12 

24 

18 

11 

25 

17 

18 

Table 6. Station-to-Satellite Range Uncertainty for Lageos. 


Based on the GEM-T2 Gravity Field Error Covariance. 




Range Uncertainty 

in millimeters 



Ascending 


Descending Overall 

Station 

Min 

Max RMS 

Min 

Max 

RMS 

RMS 

Arequipa 

7 

14 


7 

14 

10 


Easter Island 

7 

19 

12 

7 

15 

11 

11 

Goldstone 

7 

15 


7 

12 

9 


Grasse 

7 

15 

11 

7 

17 

11 

11 

Graz 

7 

15 


7 

17 

12 

11 

Greenbelt 

7 

14 


7 

14 






Table 6. Continued. 


Range Uncertainty in millimeters 
Ascending Descending Overall 


Station 

Min 

Max 

RMS 

Min 

Max 

RMS 

RMS 

Haleakala 

8 

17 

12 

8 

13 

10 

11 

Huahine 

9 

19 

13 

7 

14 

10 

12 

Kwajalein 

9 

18 

12 

8 

14 

10 

11 

Matera 

7 

14 

10 

7 

18 

12 

11 

McDonald Observatory 

7 

14 

10 

7 

12 

9 

10 

Monument Peak 

7 

15 

10 

7 

12 

9 

10 

Orroral 

7 

17 

12 

7 

15 

11 

11 

Platteville 

7 

14 

10 

7 

13 

10 

10 

Quincy 

7 

15 

11 

7 

13 

10 

10 

Greenwich Observatory 

7 

15 

11 

7 

15 

11 

11 

Shanghai 

7 

16 

11 

8 

16 

11 

11 

Simosato 

8 

16 

11 

8 

14 

11 

11 

Wettzell 

7 

15 

11 

7 

17 

12 

11 

Yaragadee 

7 

15 

11 

7 

15 

11 

11 


These results demonstrate the improvement in gravity field modeling provided by 
GEM-T2 versus GEM-T1. The overall uncertainty in the overall station-to-satellite 
range uncertainty decreased from approximately 20 mm to 10 mm. In all cases 
the range uncertainty tends to be largest when the satellite is low on the horizon 
(relative to the tracking site) and smallest when the satellite is directly overhead of 
the tracking site. This is not unexpected since the transverse component of the orbital 
perturbations are much larger than the radial component. And it is the transverse 
component that is predominantly being observed by the ranging system when the 
satellite is low in the sky. When the satellite is more directly overhead a more direct 
detection of the radial orbit error is being made. This geographic variation of the 
range uncertainty is clearly shown in the figures of the Appendix. 


8 Conclusions 

The linear perturbation approach, patterned after the method of Kaula [1966], is 
a powerful analysis technique for evaluating the accuracy of gravity field modeling. 
Given the expected uncertainty in the spherical harmonic coefficients of a gravity 
field model, via the error covariance matrix from least squares, the corresponding 
expected uncertainty in the orbital motion can be determined in a straightforward 
manner. And these results can then be extended to determine how these orbital 
modeling errors will be observed through the satellite tracking system. 

Analysis of the Lageos satellite and the gravity field models GEM-T1 and GEM- 
T2 have shown that the orbital motion is being modeled to the 1-2 cm accuracy level 
in terms of station- to-satellite range measurements. This accuracy is compatible (in 
an RMS sense) with the current accuracy of laser range measurements. The range 
uncertainty as a function of geographic location, shows that gravity field mismodeling 
will still at times be significantly above the laser range measurement accuracy, and 
thus still can be a limiting factor in Lageos geophysical studies. 
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9 Appendix 

The attached figures illustrate the variation of the range uncertainty for each of the 
tracking sites that have been considered. Figures are only presented for the case of 
using the GEM-T1 gravity field model. The results for the GEM-T2 model are very 
similar in character and primarily differ by having an overall smaller magnitude. For 
each of the tracking sites two figures are presented — one figure for the situation of 
Lageos being tracked while on an ascending pass and the corresponding case for when 
Lageos is being tracked on a descending pass. The contours of the range uncertainty 
are computed over a slightly larger geographic extent than the results provided in 
Tables 5 and 6. Here the range uncertainty is evaluated for station- to-satellite ele- 
vations down to 10 degrees (instead of 20 degrees). This larger geographic coverage 
provides for a slightly better visualization of the character of the range uncertainty. 



MINIMUM = 11 MM MAXIMUM = 26 MM RMS = 1 7 MM 
CONTOUR INTERVAL = 1 MM 




MINIMUM = 12 MM MAXIMUM = 26 MM RMS = 18 MM 
CONTOUR INTERVAL = 1 MM 





RANGE STANDARD DEVIATION PREDICTED BY GEM-T1 
EASTER IS. TRACKING LAGEOS ALONG ASCENDING PASSES 



MINIMUM = 12 MM MAXIMUM = 35 MM RMS = 21 MM 
CONTOUR INTERVAL = 1 MM 




RANGE STANDARD DEVIATION PREDICTED BY GEM-T1 
EASTER IS. TRACKING LAGEOS ALONG DESCENDING PASSES 



MINIMUM = 12 MM MAXIMUM = 31 MM RMS = 18 MM 
CONTOUR INTERVAL = 1 MM 







MINIMUM = 12 MM MAXIMUM = 26 MM RMS = 18 MM 
CONTOUR INTERVAL = 1 MM 



RANGE STANDARD DEVIATION PREDICTED BY GEM-T1 
GRASSE TRACKING LAGEOS ALONG ASCENDING PASSES 



MINIMUM = 11 MM MAXIMUM = 31 MM RMS = 19 MM 
CONTOUR INTERVAL = 1 MM 





RANGE STANDARD DEVIATION PREDICTED BY GEM-T1 
GOLDSTONE TRACKING LAGEOS ALONG DESCENDING PASSES 



MINIMUM = 11 MM MAXIMUM = 23 MM RMS = 16 MM 
CONTOUR INTERVAL = 1 MM 





RANGE STANDARD DEVIATION PREDICTED BY GEM-T1 
GRASSE TRACKING LAGEOS ALONG DESCENDING PASSES 



MINIMUM = 12 MM MAXIMUM = 37 MM RMS = 23 MM 
CONTOUR INTERVAL = 1 MM 




RANGE STANDARD DEVIATION PREDICTED BY GEM-T1 
GRAZ TRACKING LAGEOS ALONG ASCENDING PASSES 



MINIMUM = 11 MM MAXIMUM = 33 MM RMS 
CONTOUR INTERVAL = 1 MM 


19 MM 



RANGE STANDARD DEVIATION PREDICTED BY GEM-T1 
GREENBELT TRACKING LAGEOS ALONG ASCENDING PASSES 



MINIMUM = 11 MM MAXIMUM = 23 MM RMS = 17 MM 
CONTOUR INTERVAL = 1 MM 




RANGE STANDARD DEVIATION PREDICTED BY GEM-T1 
GRAZ TRACKING LAGEOS ALONG DESCENDING PASSES 



MINIMUM = 12 MM MAXIMUM = 37 MM 
CONTOUR INTERVAL = 1 MM 


RMS = 23 MM 




RANGE STANDARD DEVIATION PREDICTED BY GEM-T1 
HALEAKALA TRACKING LAGEOS ALONG ASCENDING PASSES 



MINIMUM = 13 MM MAXIMUM = 34 MM RMS = 22 MM 
CONTOUR INTERVAL = 1 MM 






MINIMUM = 11 MM MAXIMUM = 24 MM 
CONTOUR INTERVAL = 1 MM 


RMS = 17 MM 




MINIMUM = 14 MM MAXIMUM = 39 MM RMS = 24 MM 
CONTOUR INTERVAL = 1 MM 




MINIMUM = 12 MM MAXIMUM = 23 MM RMS = 16 MM 
CONTOUR INTERVAL = 1 MM 




RANGE STANDARD DEVIATION PREDICTED BY GEM-T1 
KWAJALEIN TRACKING LAGEOS ALONG ASCENDING PASSES 



MINIMUM = 15 MM MAXIMUM = 36 MM RMS = 23 MM 
CONTOUR INTERVAL = 1 MM 






MINIMUM = 13 MM MAXIMUM = 27 MM RMS = 18 MM 
CONTOUR INTERVAL = 1 MM 




MINIMUM = 12 MM MAXIMUM = 31 MM RMS = 19 MM 
CONTOUR INTERVAL = 1 MM 



MINIMUM = 12 MM MAXIMUM = 23 MM RMS = 17 MM 
CONTOUR INTERVAL = 1 MM 






RANGE STANDARD DEVIATION PREDICTED BY GEM-T1 
MATERA TRACKING LAGEOS ALONG DESCENDING PASSES 



MINIMUM = 12 MM MAXIMUM = 38 MM RMS = 23 MM 
CONTOUR INTERVAL = 1 MM 



MINIMUM = 12 MM MAXIMUM = 24 MM RMS = 17 MM 
CONTOUR INTERVAL = 1 MM 



RANGE STANDARD DEVIATION PREDICTED BY GEM-T1 
MONUMENT PK TRACKING LAGEOS ALONG ASCENDING PASSES 



MINIMUM = 12 MM MAXIMUM = 27 MM RMS 
CONTOUR INTERVAL = 1 MM 


18 MM 





RANGE STANDARD DEVIATION 
MCDONALD TRACKING LAGEOS 


PREDICTED BY GEM-T1 
ALONG DESCENDING PASSES 



MINIMUM = 12 MM MAXIMUM = 23 MM RMS = 16 MM 
CONTOUR INTERVAL = 1 MM 



MINIMUM = 12 MM MAXIMUM = 33 MM RMS = 22 MM 
CONTOUR INTERVAL = 1 MM 



RANGE STANDARD DEVIATION PREDICTED BY GEM-T1 
MONUMENT PK TRACKING LAGEOS ALONG DESCENDING PASSE 



MINIMUM = 12 MM MAXIMUM = 23 MM RMS = 16 MM 
CONTOUR INTERVAL = 1 MM 





MINIMUM = 11 MM MAXIMUM = 24 MM RMS = 18 MM 
CONTOUR INTERVAL = 1 MM 









RANGE STANDARD DEVIATION PREDICTED BY GEM-T1 
QUINCY TRACKING LAGEOS ALONG ASCENDING PASSES 



MINIMUM = 12 MM MAXIMUM 
CONTOUR INTERVAL = 1 MM 


27 MM 


RMS 


19 MM 






RANGE STANDARD DEVIATION PREDICTED BY GEM-T1 
PLATTEVILLE TRACKING LAGEOS ALONG DESCENDING PASSE 



MINIMUM = 11 MM MAXIMUM = 24 MM RMS = 17 MM 
CONTOUR INTERVAL = 1 MM 



RANGE STANDARD DEVIATION PREDICTED BY GEM-T1 
RGO TRACKING LAGEOS ALONG ASCENDING PASSES 



MINIMUM = 1 1 MM MAXIMUM = 31 MM RMS 
CONTOUR INTERVAL = 1 MM 


19 MM 






RANGE STANDARD DEVIATION PREDICTED BY GEM-T1 
QUINCY TRACKING LAGEOS ALONG DESCENDING PASSES 



MINIMUM = 11 MM MAXIMUM = 26 MM RMS = 17 MM 
CONTOUR INTERVAL = 1 MM 





MINIMUM = 12 MM 
CONTOUR INTERVAL 


MAXIMUM = 32 MM 
= 1 MM 


RMS = 21 MM 




MINIMUM = 1 1 MM MAXIMUM = 34 MM RMS = 22 MM 
CONTOUR INTERVAL = 1 MM 



RANGE STANDARD DEVIATION PREDICTED BY GEM-T1 
SIMOSATO TRACKING LAGEOS ALONG ASCENDING PASSES 



MINIMUM = 12 MM MAXIMUM = 34 MM RMS = 22 MM 
CONTOUR INTERVAL = 1 MM 



RANGE STANDARD DEVIATION PREDICTED BY GEM-T1 
SHANGAHI TRACKING LAGEOS ALONG DESCENDING PASSES 



MINIMUM = 12 MM MAXIMUM = 29 MM RMS = 20 MM 
CONTOUR INTERVAL = 1 MM 




RANGE STANDARD DEVIATION PREDICTED BY GEM-T1 
WETTZELL TRACKING LAGEOS ALONG ASCENDING PASSES 



MINIMUM = 11 MM MAXIMUM = 33 MM RMS 
CONTOUR INTERVAL = 1 MM 


19 MM 





MINIMUM = 12 MM MAXIMUM = 31 MM RMS = 19 MM 
CONTOUR INTERVAL = 1 MM 



RANGE STANDARD DEVIATION PREDICTED BY GEM-T1 
YARAGADEE TRACKING LAGEOS ALONG ASCENDING PASSES 



MINIMUM = 11 MM MAXIMUM = 26 MM RMS = 19 MM 
CONTOUR INTERVAL = 1 MM 




RANGE STANDARD DEVIATION PREDICTED BY GEM-T1 
WETTZELL TRACKING LAGEOS ALONG DESCENDING PASSES 



MINIMUM = 12 MM MAXIMUM = 37 MM 
CONTOUR INTERVAL = 1 MM 


RMS = 23 MM 



RANGE STANDARD DEVIATION PREDICTED BY GEM-T1 
YARAGADEE TRACKING LAGEOS ALONG DESCENDING PASSES 



MINIMUM = 11 MM MAXIMUM = 26 MM RMS = 18 MM 
CONTOUR INTERVAL = 1 MM 





